Realistic continuous- variable quantum teleportation with non-Gaussian resources 
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We present a comprehensive investigation of nonideal continuous-variable quantum teleportation 
implemented with entangled non-Gaussian resources. We discuss in a unified framework the main 
decoherence mechanisms, including imperfect Bell measurements and propagation of optical fields 
in lossy fibers, applying the formalism of the characteristic function. By exploiting appropriate 
displacement strategies, we compute analytically the success probability of teleportation for input 
coherent states, and two classes of non-Gaussian entangled resources: Two-mode squeezed Bell-like 
states (that include as particular cases photon-added and photon-subtracted de-Gaussified states), 
and two-mode squeezed cat-like states. We discuss the optimization procedure on the free parameters 
of the non-Gaussian resources at fixed values of the squeezing and of the experimental quantities 
determining the inefficiencies of the non-ideal protocol. It is found that non-Gaussian resources 
enhance significantly the efficiency of teleportation and are more robust against decoherence than 
the corresponding Gaussian ones. Partial information on the alphabet of input states allows further 
significant improvement in the performance of the nonideal teleportation protocol. 

PACS numbers: 03.67.Hk, 42.50.Ex, 42.50.Dv 



I. INTRODUCTION 

In recent years, theoretical and experimental investi- 
gation has been focused on the engineering of highly non- 
classical, non-Gaussian states of the radiation field (for 
a review, see e.g. 0). Interest in the production of non- 
Gaussian optical states is due to their strongly nonclas- 
sical properties, such as entanglement and negativity of 
the quasi-probability phase-space distributions, that are 
important for the efficient implementation of quantum 
information and communication protocols [H, 0, 1U, EL HIj 
and for quantum estimation tasks Q . 

Several schemes for the generation of non-Gaussian 
states, both single-mode and two-mode, have already 
been proposed @, H H EE EH, EH EH > an d many suc- 
cessful and encouragin g e xperimental realizations have 
been reported recently [H 111 EE ES El El ■ In line of 
principle, a very important result concerns the rigorous 
proof that various nonclassical properties are minimized 
by Gaussian states [20]. Therefore, it is reasonable to 
expect that the use of non-Gaussian resources may im- 
prove significantly the performance of quantum informa- 
tion protocols. In particular, concerning quantum tele- 
portation with continuous variables (CV), it has been 
shown that the success probability of teleportation can 
be greatly increased by using entangled non-Gaussian re- 
sources in the framework of the ideal Braunstein-Kimble 

(b-k) protocol i, mi, m m, m, m . 

Indeed, it has been shown that some specific two-mode 
non-Gaussian states, dubbed squeezed Bell-like states 
(that include as subcases photon-added and photon- 
subtracted de-Gaussified states (HI), when used as en- 



tangled resources provide a significant increase in the 
teleportation fidelit y o f single-mode input states under 
the ideal protocol [24(. Such an enhancement is due 
to a balancing of three different features (24| : The en- 
tanglement content of the resources, their (appropriately 
defined) degree of affinity with the two-mode squeezed 
vacuum, and their (suitably measured) amount of non- 
Gaussianity. For the precise definition of the last two 
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quantities, see Refs. 12J,|26(). It has been suggested 
that such states can be produced by combining simulta- 
neous phase-matched multiphoton processes and condi- 
tional measurements. 

The analysis of Ref. (24| has been later extended to 
consider other classes of non-Gaussian resources, such 
as two-mode squeezed symmetric superpositions of Fock 
states and of squeezed cat-like states, that allow high 
levels of performance in the teleportation of single-mode 
input states [2?| . A partial preliminary analysis of non- 
ideal cases has also been performed [271 ] by considering 
simple superpositions of independently generated fields 
converging on a common spatial volume, such as su- 
perpositions of a two-mode pure non-Gaussian resource 
and a two-mode thermal state (2^|. In this elementary 
instance, mixed non-Gaussian entangled states remain 
preferred resources for teleportation when compared to 
mixed twin-beam Gaussian states [27j |. 
In this work, using the formalism of the characteris- 
tic function, we study in full generality the Braunstein- 
Kimble protocol for C V teleportation in realistic condi- 
tions and with non-Gaussian entangled resources. Wc 
include in our investigation the main sources of deco- 
herence that lead to the degradation of the transferred 
quantum information, such as losses due to imperfect ho- 
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modyne measurements, and damping due to the propa- 
gation of the optical fields in lossy fibers. The effects of 
these inefficiencies have already been considered, among 
others, in Refs. (2^, [3(j. In particular Ref. [2^| is con- 
cerned with the study of imperfect Bell measurements, 
while in Ref. [3(3| the authors investigate the limits of 
quantum teleportation due to photon absorption during 
propagation in fibers. Besides considering each problem 
separately, these and related works are always restricted 
to the use of Gaussian resources. The main object of the 
present work is to investigate the effect of the simultane- 
ous presence of all sources of imperfection on the perfor- 
mance of CV teleportation protocols with non-Gaussian 
resources, and their robustness against decoherence. 

A general and exhaustive analysis turns out to be pos- 
sible in the framework of the characteristic function rep- 
resentation. This method has been discussed in full gen- 
erality for the description of ideal CV teleportation [3l| , 
and applied first to the case of Gaussian (3lj and non- 
Gaussian resources 0, HtJ- We will then extend the 
formalism to include the description of nonideal CV tele- 
portation, including realistic Bell measurements and de- 
coherence due to propagation in noisy channels. In order 
to investigate different optimization strategies of the non- 
ideal protocol, we will discuss optimization over the free 
parameters of the non-Gaussian resources as well as over 
the gain factor associated with the transmitted classi- 
cal information [32j (for various strategies of gain tuning 
and of optimal gain see also Refs. [3(| [H, Hil ) • Indeed, 
in the instance of non-Gaussian resources, the gain can 
be considered as a further free parameter suitable for op- 
timization. 

The paper is organized as follows. In section |TT] we 
extend the characteristic function formalism to include 
the case of realistic CV quantum teleportation. In sec- 
tion IIIII we introduce and discuss the main properties 
of some classes of non-Gaussian entangled resources. In 
section [IV] we study the efficiency of the quantum tele- 
portation protocol in the instance of fixed given values 
of the gain. In section [V] we carry out an optimization 
procedure of the protocol over the entangled resources 
and the gain parameter. Finally, in section IVII we draw 
our conclusions and discuss some outlook on current and 
future research. 



II. NONIDEAL CV TELEPORTATION 
PROTOCOL IN THE CHARACTERISTIC 
FUNCTION FORMALISM 

In this section, we describe the realistic B-K CV tele- 
portation protocol in the formalism of the characteris- 
tic function. Although several alternative formalisms are 
available for the description of the B-K CV teleportation 
protocol [H, HE H3, US HIj S3| > the characteristic func- 
tion representation proves to be particularly convenient 
when considering the nonideal case and non-Gaussian re- 
sources. The description of nonideal teleportation re- 
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FIG. 1: (Color online) Pictorial representation of the nonideal 
B-K CV quantum teleportation protocol: In the first step, 
the input mode is mixed by Alice with one of the two beams 
(modes) of the entangled resource; the ensuing state is then 
subject to a realistic Bell measurement. The result of the 
measure is communicated to Bob through a classical channel. 
In the second step, a unitary transformation, determined by 
the previous measurement, is applied to the second mode of 
the entangled resource, that is affected by decoherence during 
the propagation in a noisy channel, e.g. a lossy fiber. The 
ensuing output state is the final teleported state. 



quires the introduction of mechanisms of loss and in- 
efficiency in the main steps of the protocol. Indeed, a 
schematic description of the nonideal protocol is depicted 
in Fig. Q] 

The input state and the entangled resource states are 
assumed to be initially pure. This is not a serious limita- 
tion because one can always map the case of a nonideal 
teleportation protocol with noisy (mixed) inputs and re- 
sources to an equivalent protocol with pure inputs and 
resources but with a correspondingly larger amount of 
noise affecting the protocol. A simple illustrative exam- 
ple of this equivalence will be discussed later on in the 
present work. 

The single-mode input state (in) is mixed to mode 1 
of the entangled resource at a beam splitter. At the first 
user's (Alice) location, a Bell measurement, consisting 
in homodyne detections, is performed on the obtained 
state of mode in and mode I. In order to describe a 
nonideal measurement, one needs to model/simulate the 
inefficiencies of the photo-detectors. A realistic detector 
can be modeled by placing a fictitious beam splitter, i.e. 
a partly transmitting mirror, in front of an ideal detector 
[4l|. Based on such a prescription, a scheme describing 
a realistic Bell measurement is shown in Fig. [5] 

After a realistic Bell measurement, the result is trans- 
mitted to the receiver (Bob) through a classical channel. 
The mode 2 of the entangled resource propagates in 
a noisy channel, as a lossy fiber, to Bob's location, 
where it undergoes a unitary displacement according the 
result of the Bell measurement It is important to note 
that the generation of the entangled resource can take 
place close to the sender, and typically very far away 
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is the characteristic function of the input state, and Xres 
is the characteristic function of the entangled resource. 
By defining the quadrature operators Xj = ^( a j + a \) 

and P i = 71 (°3 - a ^ ( j = 
dated phase-space real variables 



in, 1,2), and the asso- 
the characteristic function can 



be written in terms of Xj, Pj, i.e. Xo{ a inl ot\\ a 2 ) = 

Xo( x in,Pm,Xx,pi;X2,P2)- 

The first step of the protocol consists in the Bell measure- 
ment at Alice's location, that is the homodync measure- 
ments of the first quadrature of the mode 1 and of the 
D second quadrature of the mode in, with results x and 

p, respectively. After such nonidcal Bell measurement, 
the remaining mode 2 is left in a mixed state described 
by the corresponding single-mode characteristic function 
XBm{x2,P2)- In appendix [Al we prove in detail how to 
compute it in full generality for arbitrary single-mode in- 
puts and arbitrary two-mode entangled resources. Here 
we report only the final expression: 



FIG. 2: (Color online) Model scheme of a realistic Bell mea- 
surement. The model takes into account the non-unity effi- 
ciency of the detectors D performing the homodyne measure- 
ment. In the depicted scheme, such inefficiency is simulated 
by the introduction of two fictitious beam splitters BS2 and 
BS3, with equal transmissivity n. 

from the receiver, as one of the main tasks of quantum 
teleportation is the transfer of quantum information 
across long distances. It is then legitimate to assume 
that the radiation field associated to mode 1 is not 
affected by losses due to propagation, while the field 
associated to mode 2, which usually has to propagate 
over much longer distances, can be strongly affected by 
decoherence. The degradation of quantum information 
is caused by the propagation of field mode 2 in a noisy 
channel. Therefore, the output teleported state depends 
both on the inefficiency of the homodyne detectors and 
the decoherence rate of the noisy channel. 

We can formalize the effects of the above-described dy- 
namics as follows. Let us denote by pi„ = \4>)inin{4'\ 
and p res = |V>)i2 12(^1 the density matrices associated, 
respectively, with the single-mode pure input state and 
with the two-mode pure entangled resource. The single- 
mode input state is initially disentangled from the two- 
mode entangled resource, so that the the initial three- 
mode field is po — Pin ® Pres , and the initial global char- 
acteristic function reads 

Xoi^in, an a 2 ) = Tr[p D in (a m ) Di(ai) D 2 (a 2 )} 

= Xin{a in ) Xres{ai;a 2 ) , (1) 

where Tr denotes the trace operation, Dj (ctj ) denotes the 
displacement operator for the mode j (j = in, 1,2), Xin 



X Bm (x 2 ,P2) = - ' ^ 2 X) J dtdv e^-™ x 

expj-^ 2 + W 2 )}, (2) 

where T 2 and R 2 = (1 — T 2 ) denote, respectively, the 
transmissivity and reflectivity of the beam splitters that 
model the losses. The function V(p, x) is the distribution 
of the measurement outcomes p and x (see appendix [A"]l . 
Note that the Gaussian exponential in Eq. |J2J) is related 
to the vacua entering the input ports of the fictitious 
beam splitters. 

Afterwards, mode 2 propagates in a damping channel, 
like, e.g., a lossy fiber, before it reaches Bob's location. 
The Markovian dynamics of a system subject to damping 
is described, in the interaction picture, by the following 
master equation for the density operator p [iH, 51] : 



dtp = — \ n t hL[a\]p+ {nth + ^)L[a 2 ]pj , (3) 

where the Lindblad supcropcrators are defined as 
L[0]p = 2e>p£> t -O t Op-pO t C, T is the mode damping 
rate, and nth is the number of thermal photons. Finally, 
at Bob's location, a displacement A — g(i + ip) is per- 
formed on mode 2. The real parameter g is the so-called 
gain factor The combined effect of propagation in 

a damping channel and unitary displacement determines 
the characteristic function Xout(x 2 ,p 2 ) of the final out- 
put state of the teleportation protocol (see appendix [5] 
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for details): 

Xout(x 2 ,p 2 ) = Xtn(gTx 2 ,gTp 2 ) x 

Xres (gTx 2 ,-gTp 2 ;e~^x 2 ,e~^p 2 ) x 

exp|-^r TiJi (a;l+^)|, (4) 

where t = Tt, and the thermal " renormalized" phase- 
space covariance T r ^ is defined as: 

Tr,B = (1 - e" T ) Q + n th \ + g 2 R 2 ■ (5) 

The form of Eq. (j4|) highlights the different roles played 
by the two sources of noise introduced in the tclcpor- 
tation protocol, associated, respectively, to the damp- 
ing rate T and the reflectivity R 2 . The two decoher- 
ence mechanisms act separately but also in combination, 
as one can see from the Gaussian exponential factor in 
Eq. (|3]), which in fact is nonvanishing for R ^ and/or 
t 7^ 0. The effect of the imperfect Bell measurement 
is expressed also by the presence of the scale factor T 
in the arguments of the input and resource characteris- 
tic functions Xin and Xres- Viceversa, decoherence due 
to noisy propagation affects obviously only mode 2 by 
means of the exponentially decreasing weight e _ 2 in the 
arguments of Xres- The factorized form of the output 
characteristic function, holding for the ideal protocol [3l| , 

Xout(x 2i p 2 ) = Xin(x 2 ,P 2 ) Xres(x 2 ,-P 2 ;x 2 ,p 2 ), (6) 

is recovered, as expected, from Eq. ([4]) when R = (T = 
1), T = (r = 0), and g = 1. 

III. ENTANGLED RESOURCES: TWO CLASSES 
OF OPTIMIZED NON-GAUSSIAN STATES 

Given the general description of the nonideal protocol 
in terms of the characteristic functions, in this section we 
analyze the performance of two classes of non-Gaussian 
entangled resources for the teleportation of input coher- 
ent states, respectively, the two-mode squeezed Bell-like 
states \ip)sB and the two- mode squeezed cat-like states 
\ip)sc- 

Wsb = S 12 (C){cos<5|0,0) + e* e sin<5|l,l)}, (7) 
\i>) sc = A/" 5C S 12 (C){cos <5|0,0) +e w sin<5| 7 , 7 }}, 

(8) 

where 5i2(C) = e~'» a i a 2 + ^ aia2 is the two-mode squeezing 
operator, £ = re 1 *, \m,n) = \m)\ <g> \n) 2 is a two-mode 
Fock state (of modes 1 and 2), 17,7) = | 7 )i <g) 17)2 is 
a symmetric two-mode coherent state with complex am- 
plitude 7 = |7|e ll/3 , and the normalization factor Msc is 



Nsc = {1 + e I tI 2 sin 28 cos #} -1 / 2 . In order to obtain a 
maximization of the teleportation fidelity, it is necessary 
to perform a simultaneous balanced optimization, on the 
free parameters (5, 6, 7), of some partially competing 
properties 0, [27| . These include the entanglement con- 
tent, the amount of non-Gaussianity of the state [2(|, and 
a squeezed- vacuum- affinity Q. For a generic pure state 
\ip) : the latter is defined as [H, H3| : 

g = sup|(-r#)| 2 , (9) 

r 

with I — r) = Si 2 (— r)|0, 0). Indeed, the optimal 
non-Gaussian resources ([7]) and (JSJ) exhibit a sufficient 
squeezed-vacuum-affinity, which then appears to be a 
crucially needed property in order to select efficient and 
highly performing non-Gaussian resources. For instance, 
one could as well consider a different form of squeezed 
Bell-like state, the so-called "Buridan donkey" or "Ham- 
let" state \ip)sB 2 , that is obtained from the singlet Bell 
state as follows: 

U') SB2 = S 12 ((){cos5\0, 1) + e w sin<S|l,0)}. (10) 

It is indeed simple to verify that, although such a state, at 
fixed squeezing, is more entangled than a Gaussian twin 
beam, it performs less efficiently both in the ideal and 
in the realistic teleportation protocol. This fact can be 
understood if one looks at the behavior of the squeezed 
vacuum affinity [24[ . Namely, the Buridan donkey state 
Eq. pop does not contain the fundamental Gaussian con- 
tribution coming from the squeezed vacuum. Therefore, 
it is "unbalanced", in the sense that it is less affine 
to the squeezed vacuum, and excessively non-Gaussian 
compared to the optimized Bell-like state \ip)sB- There- 
fore, the fine interplay among these three quantities, i.e. 
the entanglement, the degree of non-Gaussianity, and, 
in particular, the squeezed vacuum affinity, cannot be 
realized in the non-Gaussian resource (TU|) [24| . The 
crucial role played by the squeezed vacuum affinity for 
the performance of different non-Gaussian resources has 
been studied in Ref. (27[. In particular, it has been 
shown that, in the ideal teleportation protocol, the two- 
mode squeezed symmetric superposition of Fock states, 
i.e. Si 2 (C) X)fc=o c k\k, k), when optimized for the telepor- 
tation of both input coherent states and single-photon 
states, reduces to a squeezed truncated twin beam, i.e. 
<Si2(— r) Y^k=o tanh fc s\k, k). In the same paper, it has 
been also shown that, in the ideal teleportation proto- 
col, the optimized two-mode squeezed cat-like states, i.e. 
Eq. JS} , possess a high amount of squeezed vacuum affin- 
ity. Therefore, besides a certain amount of entanglement, 
also a sufficient degree of squeezed vacuum affinity ap- 
pears to be necessary for a non-Gaussian resource to be 
optimal for a B-K teleportation protocol. 

Clearly, the performance of B-K teleportation proto- 
cols depends strongly on the structure of the second- 
order correlations in the entangled resources. In this 
sense the B-K protocol, with its structure of homodync 
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measurements, is particularly tailored to the use of Gaus- 
sian resources. Therefore, a non-Gaussian resource may 
improve on the performance of a corresponding Gaus- 
sian one only if the fundamental Gaussian contribution 
coming from the squeezed vacuum is subject to a not 
too drastic modification. A large value of the affinity 
assures that the non-Gaussian resource satisfies such a 
requirement. The interplay of the affinity with the non- 
Gaussianity and the degree of entanglement allows to sin- 
gle out those non-Gaussian resources possessing higher- 
order correlations that add to the leading Gaussian struc- 
ture of the two-mode entangled resource, thus enhancing 
further the protocol efficiency, and lacking those non- 
Gaussian contributions that are incompatible with the 
structure of the B-K protocol. 

Indeed, a very interesting question open for future in- 
vestigation is the inverse of the one studied in the present 
paper. Here we are analyzing the problem of optimizing 
non-Gaussian resources given the B-K protocol. The in- 
verse question would be that of adapting the protocol 
to the resources. Namely, given a certain class of non- 
Gaussian squeezed resources with some given properties, 
one asks how the B-K protocol would have to be modified 
in order to optimize the fidelity of telcportation. 

We now proceed to determine the general expression 
of the fidelity of teleportation in terms of the character- 
istic function for the three different non-Gaussian entan- 
gled resources \tp)sB, W)sc, \4'}sb2- For instance, the 
two- mode characteristic function xsb of the entangled 
resource (0 reads 

Xsb{oci, a 2 ) = Tr[\ip)sB sb(M Di(ai)D 2 (a 2 )} , (11) 

and analogous expressions hold for xsc and XSB2- The 
corresponding explicit expression is obtained using the 
two-mode Bogoliubov transformations 

^12 (0 a i Si2(C) = coshr — e 1 ^ sinh raj-, 

(Mi = 1,2), (12) 

and the relation 

(m\D(a)\n) = ' ' a m ~ "e^L^ n \\a\ 2 ) , 

(13) 

where L„ (•) denotes the associated Laguerre poly- 
nomial of order n. 

The quantity measuring the success probability of a 
teleportation protocol is the fidelity of teleportation T = 
Tr[pi n p out ]. In the formalism of the characteristic func- 
tion, the fidelity reads 

F = — J d 2 a Xin(a)xout(-a) , 

— J dx 2 dp 2 Xin(x2,P2)Xout(~X 2 , -p 2 ) , (14) 

where a = -^(x 2 + ip 2 ), d 2 a = \dx 2 dp 2 , and Xout{oi) = 
Xout (x 2 ,p 2 ) is given by Eq. (0}. In the case of input co- 
herent states pi n = \{3}i n in{P\ with complex amplitude 



f3, that we will always consider in the following, the char- 
acteristic function of the input Xin( a ) reads: 

Xin (a) = e -§l«l 2 +M*-^) . (15) 

Eq. (|14[) is the fundamental quantity that measures the 
efficiency of a CV teleportation protocol (ideal or non- 
ideal). At fixed squeezing, the optimization procedure 
consists in the maximization of the teleportation fidelity 
(TTJJ over the free parameters of the non-Gaussian entan- 
gled resources, Eqs. ([7]), ©, and (TIT))) . Since it can be ver- 
ified explicitly that the optimal choice for the phases <j>, 
0, and (p are <j> = ir and 6 = ip = 0, the squeezed Bell-like 
state \iPsb) and the Buridan donkey state \^SB2) have 
a unique available free parameter S. On the other hand, 
the squeezed cat-like state \ipsc) has two free parame- 
ters, the angle <5 and the modulus |7|. The analytical ex- 
pressions of the telcportation fidelities of input coherent 
states corresponding to the three different classes of non- 
Gaussian resources, respectively F^ir, •^sc( r ' ^> M)' 
and J~g 9 g2{r,S) are reported in appendix [Cl Eqs. (|C1[) . 
(|C4[) . and (|C6[) . Let us notice that we have introduced 
the superscript (<?) to explicitly indicate the dependence 
on the gain g. It is simple to verify that, for arbitrary g, 
the fidelities are explicitly dependent on the amplitude 
(3 of the input coherent states. In the next sections, the 
(numerical) optimization procedures of the fidelity T will 
be implemented following two different routes. In section 
|TV]we operate at a specific value of the gain g = 1/T, 
for a fixed value of the transmissivity T 2 . This is the 
only choice that makes the fidelity independent of (3. In 
section [V] wc adopt a more general approach by letting 
j be a fully free parameter and performing appropriate 
optimization procedures. 

In both cases, the maximization is carried out at fixed 
(finite) squeezing r, and at fixed r, n t h, and R. From 
an operational point of view, fixing these parameters is 
equivalent to assume control on the characteristics of the 
experimental apparatus, including the inefficiency of the 
photo-detectors and the length and damping rate of the 
noisy channel. Finally, concerning the experimental real- 
ization of the two-mode non-Gaussian resources, Eqs. ([7]) 
and (151, a detailed theoretical proposal is put forward in 
Ref. [241 ] . The experimental realization of the single- mode 
version of the state \ip)sc has been reported in Ref. [lj|. 

IV. /3-INDEPENDENT OPTIMAL FIDELITY 

In this section, we analyze the success probability of 
quantum teleportation for the gain g fixed to be g = 
1/T. It is immediate to verify that with this choice, the 
fidelity becomes /3- independent (see appendix ICj) . This 
choice allows to assume no knowledge about the alpha- 
bet of input coherent states, while in the next section 
we will assume a partial knowledge on the input states 
over an interval of values of (3. For g = 1/T, the ex- 
pressions for the fidelities ^sb^jS), Tsc ( r i S, and 
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3~SB2(r, 8) (where the superscript (g) has been removed) 
greatly simplify. At fixed r, r, nth, and R, the optimal 
fidelities of teleportation are defined as: 

F^t = max J SB (r,(J), (16) 

fip? 2) = max FsB2(r, 5), (17) 

^ t C) = mBxT sc (r,6,\7\). (18) 
s,h\ 

In Fig. [3]we plot the optimal fidelities, corresponding to 
the three classes of non-Gaussian resources, as functions 
of the squeezing r at different values of the parameters r, 
nth, and R = \fl —T 2 . For comparison, the fidelity asso- 
ciated with the Gaussian squeezed vacuum (twin beam) 
is reported as well. In order to understand the separate 
effects of the two different sources of decoherence on the 
degradation of the fidelity, we consider two cases: (i) de- 
coherence due to imperfect Bell measurements alone, i.e. 
R > and t = (see Fig. [3] panel I); (ii) decoherence 
due to propagation in noisy channels alone, i.e. R = 
and r > (see Fig.[3]panel II). 

In the first case, the fidelities grow monotonically, with 
increasing r, tending towards an asymptotic saturation 
value. This behavior is equivalent to that observed in 
the instance of an ideal protocol with noisy resources 
[27j . Indeed, the case of a nonideal teleportation protocol 
with noisy (mixed) resources is equivalent to the case of 
a nonideal protocol with pure resources but with a larger 
amount of noise. 

In the second case, as r increases, the fidelity first in- 
creases up to a T-dependent maximum r max (r) , and then 
decreases for larger values of r. This behavior can be ex- 
plained observing that there are two competing effects 
associated to increasing the degree of squeezing. The 
first effect is constructive and is due to the enhanced 
affinity of the entangled resource with an EPR state for 
increasing r. This constructive effect is contrasted by a 
disruptive one due to the optical photons generated by 
the squeezing that add to the thermal photons of the 
channel (initially set to zero). For not too large values of 
r, the first effect dominates, until a maximum is reached 
at r = r max (j). For r > r max {r) the disruptive effect 
becomes dominant, and the increasingly large number of 
optical photons amplifies the decoherence, leading to a 
strong suppression of the fidelity. The interplay between 
squeezing r and channel decay rate r can be understood 
quantitatively by investigating the structure of the out- 
put characteristic function Xout Eq. (U]) that enters in the 
expression of the fidelity (fl4|> . For gT = 1, Xout takes the 
form 

Xout{x 2 ,P2) = e-' r - B ^ + ^ ) X m (x 2 ,p 2 ) X 

where T Ti e is given by Eq. ([5]). We see that, if r ^ 0, 
the exponential weights e~ T /' 2 introduce an asymmetry 



between the two modes of the resource in the expression 
of the characteristic function Xres ■ This asymmetry is re- 
sponsible for the decrease of the fidelity for r > r max {T) 
at t ^ 0. The important ensuing conclusion is that it 
is in fact detrimental to increase too much the squeezing 
when the losses cannot be reduced strongly. Therefore, 
the primary experimental goal should always be that of 
reducing the losses rather than incrementing the squeez- 
ing. 
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FIG. 3: (Color online) Optimal fidelities of teleportation J- op t, 
as functions of the squeezing parameter r, for different val- 
ues of the parameters r, nth, and R. The fidelities cor- 
respond to the teleportation of single-mode input coherent 
states |/3) using two-mode squeezed Bell-like states (full line) 
or two-mode squeezed cat-like states (dashed line) as entan- 
gled resources. The fidelities associated with the two-mode 
squeezed vacuum (dotted line) and to the Buridan donkey 
states (long-dotted line) are reported for comparison. In panel 
I, r = 0, nth = 0, and the reflectivity is fixed at the values 
R 2 — 0, 0.05, 0.1, 0.15. For each entangled resource (associ- 
ated with a specific plot style), the corresponding curves are 
ordered from top to bottom with increasing R 2 . In panel II, 
nth = 0, R = 0, and the reduced time is fixed at the values 
r = 0, 0.1, 0.2, 0.3. For each entangled resource (associated 
to a specific plot style) the corresponding curves are ordered 
from top to bottom with increasing r. 

Moreover, Fig. [5] (Panel II) shows that indeed at cur- 
rent experimentally attainable values of the squeezing, 
i.e. r < 1.5, the nonideal teleportation protocol operates 
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already in the regime of best efficiency, and both the 
squeezed Bell-like resources J7]) and the squeezed cat-like 
resources (JHJ) perform much better than the correspond- 
ing (i.e. at the same squeezing) Gaussian resources. This 
result generalizes and confirms the analogous behavior 
observed in the instance of ideal protocols [HI, H3 ■ On 
the contrary, as already anticipated in the previous sec- 
tion, the Buridan donkey resources allow for teleporta- 
tion fidelities even worse than those associated with the 
Gaussian twin beam. Indeed, from Eq. (|CJ6|1 it follows 



that the optimal /3-independent fidelity F^pf 2 ^ is ob- 
tained (with g = 1/T) by letting (5 = 0. In this case, 
the Buridan donkey state trivially reduces to a two-mode 
squeezed Fock state. 

It is worth noting that, at gT = 1 and any fixed value of 
r, all the non-Gaussian resources share with the Gaussian 
one the same maximum value of the fidelity, obtained at 
the same value r max (r) of the squeezing parameter. This 
implies that, at given r, one can determine the value 
r mnl (r) for all the various resources by just considering 
the simple Gaussian instance. A straightforward compu- 
tation then yields 



cxp{2 T raax } 



/cosh(r/2) + 1 
cosh(r/2) - 1 



(19) 



We see that, for increasing r, the range [0, r max (r)] of 
best efficiency reduces. 




FIG. 4: (Color online) Optimal fidelities of teleportation T ov t, 
as functions of the squeezing parameter r, with r = 0.3, 
n t h = 0, and R 2 = 0.05. The different lines represent the 
fidelity of teleportation of input coherent states ]/3), corre- 
sponding, respectively, to the following entangled resources: 
Squeezed Bell- like state(full line), squeezed cat-like state 
(dashed line), squeezed vacuum (dotted line), Buridan don- 
key state (long-dotted line), and photon-subtracted squeezed 
state (dot-dashed line). 

Finally, we consider the combined effect of the two de- 
coherence mechanisms. In Fig. [4] we plot the optimal fi- 
delities associated to the different classes of non-Gaussian 
resources with the experimental parameters fixed at the 
values t = 0.3, nth — 0, and R 2 = 0.05. In this case, the 



simultaneous presence of the two effects leads to a strong 
suppression of the fidelity both with respect to the ideal 
case and to each of the two nonideal cases taken sepa- 
rately. A regime of best efficiency is still present, but 
significantly reduced. 

In Fig.[3] we also report the fidelity of teleportation as- 
sociated with the two-mode photon-subtracted squeezed 
states: 



IPSS 



AAa ia2 S 12 (C)|0,0) 



= A^Si 2 (C){-|0,0)+e^tanhr|l,l)} , (20) 

where N = (1 + tanh 2 r) -1 / 2 is the normalization [24| . 
This non-Gaussian resource belongs, as a particular sub- 
case, to the class of the squeezed Bell-like resources. In- 
deed, Eq. ([7]) reduces to Eq. (j2"0)) for 5 = arctan(tanh r) 
(with = 7T and 6 — 0). The interest in the resources 
([2T))) is due to the fact that such states have been already 
produced in the laboratory 0, [T^j- Furthermore, the 
corresponding de-Gaussification scheme can be easily in- 
tegrated in the B-K teleportation protocol [2l|. As in the 
ideal instance [24]], also in the realistic case, see Fig. 0] 
the performance of the non-Gaussian resource (|20p is in- 
termediate between that of the Gaussian twin beam and 
of the squeezed Bell-like states for < r < 1, while, 
for r > 1, it degrades much faster. The affinity to the 
two-mode squeezed vacuum of state ([213) decreases for 
growing r; correspondingly, the resource becomes more 
and more non-Gaussian. Moreover, for some intervals of 
values of the squeezing parameter the photon-subtracted 
squeezed states behave better than the squeezed cat-like 
states. It is worth noticing that there exists a "crossing" 
value of r at which the optimal squeezed Bell-like state 
reduces to the photon-subtracted squeezed state and thus 
the corresponding fidelities of teleportation coincide, see 

Fig.m 

In this section we have considered always the case 
gT = 1. The scenario changes dramatically if gT ^ 1. 
Indeed, in this case the analytical expressions of the fi- 
delities JC1|), ([C4|), |C5|, Jc||, depend on the coherent 
amplitude [3. This dependence affects quite significantly 
the behavior of the fidelities as we will see in the next 
section. 



V. AVERAGE OPTIMAL FIDELITY AND 
ONE-SHOT FIDELITY 

In order to investigate possible improvements in the 
efficiency of the teleportation protocol, in this section we 
aim at optimizing the success probability of teleportation 
assuming g as a further free optimization parameter. The 
fidelity of teleportation as a function of the gain is studied 
for several input states in Ref. 
strategies are considered in Ref 



33 , while displacement 
34j | , in order to improve 



the output quality for a reduced alphabet of possible in- 
put states. These two important works are concerned 
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with the study of the ideal protocol implemented using 
Gaussian resources. The effect of absorption due to prop- 
agation in fibers is studied in Ref. [3(| , where, for the case 
of Gaussian resources, it is shown that the gain-optimized 
fidelity of teleportation is strongly suppressed. 
Let us now describe the optimization procedure applied 
to the instance of non-Gaussian resources. Following the 
approach of Refs. [10, [H[ , we define the average fidelity 
T by averaging the /3-dcpendcnt fidelity T{f3) over the 
set of input coherent states \(3) as follows: 



T 



d 2 {3T{{3)p{{3) 



(21) 



For each possible value of (3 one can estimate the suc- 
cess probability of teleportation associated with a spe- 
cific event. The functions J-i s (f3,r) yield the telepor- 
tation fidelities at given squeezing r and for an input 
coherent state with specific amplitude (3. Partial infor- 
mation about the alphabet of input states, quantified by 
the choice of the variance a in the distribution (|22|) . can 
be exploited to obtain a refinement of the optimization 
procedure. Indeed, we expect that smaller values of a, 
corresponding to a better knowledge of the alphabet of 
input states, will lead to higher values of the one-shot 
fidelities. 



p{(3) = (na) 1 cxp{— er 



(22) 



where the function p{(3) (|2"2"|) is a Gaussian distribution 
centered at (3 = 0. The variance parameter a determines 
the cutoff of the amplitude /3, and thus the reduced alpha- 
bet that one considers. We will compare our results with 
the quantum benchmark for the storage and transmission 
of coherent states distributed according to Eq. (|2"2"|) [44| . 
This benchmark is equivalent to the upper bound T c iass 
achievable with any classical strategy, and satisfying the 
inequality j44|: 



•J i clas 



< 



1 



2(7+1 



(23) 



The first step of our optimization procedure, employing 
Eqs. (jCip and (|C4[) . is to determine the average fidelities 



^sb{t, S) and ^^(r, 5, 17Q, whose expressions we do not 
report because their long and cumbersome structure is 
not particularly illuminating. Further, we do not apply 
the optimization procedure to the teleportation fidelity 
associated to Buridan donkey resource as we have already 
showed that no enhancement can be obtained compared 
to the Gaussian twin beam resource. At fixed squeezing r 
and experimental parameters r, n t h, and R, we define the 
optimal values <? op t, S op t, and |7 op t| of the free parameters 
as those that maximize the fidelities averaged over the 
values of f3 weighted according to the normal distribution 



FsT \ r ^o P t) = maxJ"/ B (r,(7), 



(24) 



Fsc Pt '(r, S opt , \lo P t\) = r max Fsc(r, 5, (25) 

The optimal values of the parameters are determined nu- 
merically. Next, we introduce the one-shot fidelities T\ s 
as the non-averaged fidelities evaluated at the optimal 
values of the parameters and a fixed value of (3: 



(26) 



^ ( f c) (/3,r) = ^ 9 S pt) (l3,r, S opt , \ 7opt \) . (27) 



A. Fidelities: Variable r, fixed r 

In Fig. El at the same fixed parameters of Fig. |4j 
t = 0.3, n t h ~ 0, R 2 = 0.05, we plot the one-shot fi- 
delity J-\ s as a function of r, both for the non-Gaussian 
resources and for the Gaussian twin beam. In panel I 
we plot the one-shot fidelities for a — 10 and (3 = 1,2,3; 
in panel II we plot the one-shot fidelities for a = 100 
and (3 = 3,5,10 (obviously \(3\ 2 must fall in the inter- 
val [0,cr]). Let us notice that according to Eq. (|23[) . 
the quantum benchmarks for the two choices of a are: 
Fcla SS (v = 10) « 0.523 and F class (a = 100) « 0.502. 

Panel I of Fig. corresponding to a = 10 and small 
values of \/3\, shows that, as soon as r is different from 
zero, all the resources yield fidelities above the quan- 
tum benchmark. Moreover, one observes a significant 
enhancement of the fidelities with respect to the (3- 
independent ones of Fig. [4] (corresponding to g = 1/T). 
Indeed, while the /3-indcpendcnt fidelities in Fig. |4] are 
well below the value 0.8, all the one-shot fidelities of panel 
I of Fig. [5] lie above this value for squeezing r ranging 
from about 0.8 to about 1.2, depending on the resource 
being considered. Also in this case the non-Gaussian re- 
sources exhibit better performances with respect to the 
Gaussian ones. Panel II of Fig. [5] shows that for a = 100 
and larger values of |/3|, the enhancement of the fidelity 
is quite modest compared to the /3-indcpcndcnt fidelity, 
reported in Fig. 21 This result is not surprising because a 
variance a = 100 obviously allows less knowledge on the 
alphabet of input states with respect to the case a = 10. 
It is important to remark that the curves corresponding 
to the same entangled resource but different values of (3 
become effectively distinguishable only r > 1. 



B. Fidelities: Variable r, fixed r 

We now study the fidelities as functions of the reduced 
time, or effective length of the fiber r, at fixed squeezing 
r. To this aim, fixing the squeezing parameter at the 
intermediate value r = 0.8, with nth = 0, R 2 = 0.05, 
we investigate the behavior of the one-shot fidelities. In 
Fig. [6j choosing the same values of (3 and a as in Fig. 
we plot J- is as a function of r. 
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FIG. 5: (Color online) One-shot fidelity of teleportation Ti s , 
for input coherent states \ j3), as a function of the squeezing pa- 
rameter r, with r = 0.3, nth = 0, R 2 = 0.05, for the following 
resources: Squeezed Bell-like state (full line); squeezed cat- 
like state (dashed line); and a squeezed vacuum (dotted line). 
In panel I /3 = 1, 2, 3 and a — 10. In panel II j3 = 3, 5, 10 
and a = 100. For each entangled resource (associated with a 
specific plot style) , the corresponding curves are ordered from 
top to bottom with increasing f3. 



FIG. 6: (Color online) One-shot fidelity of teleportation T\ a , 
for input coherent states \f3), as a function of the reduced 
time r, with r = 0.8, n t h = 0, R 2 = 0.05, for the follow- 
ing entangled resources: Squeezed Bell-like state (full line); 
squeezed cat-like state (dashed line); and squeezed vacuum 
(dotted line). In panel I = 1, 2, 3 and a — 10. In panel 
II /3 = 3, 5, 10 and a — 100. For each entangled resource 
(associated to a specific plot style), the corresponding curves 
are ordered from top to bottom with increasing /3. 



VI. CONCLUSIONS 



Fig. [6] shows that the teleportation fidelity remains 
above the classical threshold up to significantly large val- 
ues of t. At a = 10 and for small values of of the coherent 
amplitude (3 (see panel I), the one-shot fidelities associ- 
ated to the same resource but to different values of j3 
arc distinguishable. Viceversa, at a = 100 and for larger 
values of (3 (see panel II), the fidelities associated to the 
same resource but different values of (3 are virtually indis- 
tinguishable, and thus effectively j3- independent. Com- 
paring the performance of the different resources in the 
time domain, one can distinguish three regimes: For 
short and long times the one-shot fidelities associated 
to non-Gaussian resources exhibit a significant enhance- 
ment compared to the Gaussian instance, while for in- 
termediate times there is a substantial equivalence in the 
performance of Gaussian and non-Gaussian resources. 



In this paper we have investigated the performance of 
non-Gaussian entangled resources in nonideal protocols 
of quantum teleportation of input coherent states. We 
have resorted to the characteristic function formalism for 
the description of protocols affected by decoherence. We 
have discussed how the effects of decoherence stemming 
from photon losses in fiber and imperfect Bell measure- 
ments affect the success probability of teleportation. In 
particular, we have established that, while the fidelities 
associated to different resources remain above the clas- 
sical benchmark for quite long times, the non-Gaussian 
resources perform always better than the Gaussian ones, 
in the ideal as well as in the nonideal quantum telepor- 
tation protocol. 

The present analysis should be extended to include 
teleportation of two-mode states using multimode non- 
Gaussian resources [45j . It would be also interesting to 



10 



consider the optimization of teleportation protocols with 
non-Gaussian entangled resources with respect to local 
properties (such as single- mode squeezing), by extending 
the existing schemes for the local optimization of Gaus- 
sian resources 14a. 
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APPENDIX A: REALISTIC TELEPORTATION 
PROTOCOL IN THE CHARACTERISTIC 
FUNCTION REPRESENTATION 

Here we outline in some details the description of 
the nonideal Braunstein-Kimblc teleportation protocol 
in terms of the characteristic function formalism. Let 
us consider the Bell measurement Fig. O The 50/50 
beam splitter BS\ acts on the initial density operator 
Po through a SU(2) transformation UbSi(j), consisting 
in a balanced mixing of the modes in and 1: 



Po = Ubs, (J) Pa Ul Sl (J 



(Al) 



This transformation yields the following relation among 
the phase-space variables corresponding to the input and 
output modes of the beam splitter: 



73 ( x * 



■x 1 



Pi 



71 fr* 



Pi) 



Xl) , Pi = -^(Pin ~ Pi) , 



(A2) 



where the canonical variables x'j and p'- correspond to the 



primed modes j' (j' = in' , 1'). Using Eq. (|A2[) . the char- 
acteristic function xb associated to the density operator 
p' reads: 

x'o{x' in ,p' %n ]x' 1 ,p' 1 - 1 x 2 ,P2) = 



Xk(x'k,P'k) 



3,4. 



(A4) 



The beam splitters BS2 and BS3, associated to the 
SU(2) transformations Ubs^t]) (j = 2, 3), transform the 
variables x'j and p'j (j = in, 1,3,4) according to the fol- 
lowing relations: 



x'{ = Tx\ - Rx' 4 , p'{ = Tp\ - Rp' A , 
x'i = Tx' 4 + Rx[ , p'l = Tp' A + Rp[ , 



(A5) 



Tx' in - Rx' 3 , p" in = Tp' in - Rp' 3 



3 ■ 



x'i = Tx' 3 + Rx' ln , p'l = Tp' 3 + Rp'„ 



(A6) 



where T 2 = 77 and R 2 = 1 — 77 denote, respectively, 
transmissivity and reflectivity of the beam splitters, and 
double primes denote output modes. At the output of 
BS2 and BS3, the transformed five-mode density opera- 
tor reads 



Pes = Ubs 3 {v)Ubs 2 {v)Pq ® l°>3' 3' (0| <8> 

\0)4'4'(0\ul S2 (n)ul S3 (n), (A7) 



( rjSf 04, + 4) , ^(4 + pi) 

1 -^( X tn ~ X 'l) > -miP'in -P'l); x 2 ,P2 ) -(A3) 



V2 X 



V2 X 



where p' is given by Eq. (|A1[) . By exploiting the above 
transformations, a simple formal manipulation yields the 
characteristic function \bs corresponding to pbs'- 
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XBS {x'L , Pi„ 5 X" ,P";X 2 ,P2\X3, p% ] x'l ,p'l) = X 



±= [Tx'[ n + Rx'i + Tx'l + Rx'l] , [Tp'[ n + Rp'i + Tp'{ + Rp'^ 



±= [Tx'{ n + - Tx'{ - Rx'l] , i= [Tpl + Rp'i - Tp'{ - Rp'i] ; x 2 , p 2 ] X 3 (Tt% - Rx'( n , Tp'i - Rp 



V2 



Xa 



(Tx'l-Rx'l,Tp'i-Rp'A , 



(A8) 



r 



where the functions Xk, with k = 3,4, are given by 
Eq. (|A4|) . The Bell measurement, that consist in the 
two homodyne measurements of the variables p" in and 
x'{, yields the results p'( n = p and x'[ = x, and trans- 
forms the characteristic function xbs to (see appendix [B] 
for details): 

XBm(x2,P2) = V ^: 2 X) f dx'l n dp'{ j*'L*-™v" X 



(27T) 2 



XBs(xi n ,0;0,p'{; x 2 ,p 2 ; 0,0; 0,0) . 



(A9) 



Here V(p, x) is the distribution function of the outcomes 
p and x. It reads: 

V(p,x) = Tr[\p) in » in "(p\ (g> \x)i" i»(x\ pbs] 



(27T) 2 



dx'LdPi exp{ix'- n p - ixp"} x 



X bs(x'L,0; 0,^:0,0; 0,0; 0,0) 



(A10) 



Let us notice that in Eq. (|A9[) the modes 3" and 4", asso- 
ciated with the unused output ports of the two fictitious 
beam splitters, have been traced out. This is equivalent 
to setting to zero in the function xbs the variables x" 
and p'- with j = 3, 4. Meanwhile, mode 2 propagates in a 
noisy channel towards Bob's location. The dynamics of 
this mode is described by the master equation ([3]), and 
the corresponding diffusion equation for the characteris- 
tic function Xt(x 2 ,p 2 ) is (42l. [43l| : 



dt Xt{x 2 ,p 2 ) 



IV 1 



+ n t h)(xl +pl) 



+x 2 d X2 +y 2 d V2 Xt[x 2 ,p 2 ) 



(All) 



Given the initial characteristic function XBm(x 2 ,p 2 ) 
Eq. (|A9|) . the global characteristic function Xt{x 2 ,p 2 ) fol- 
lowing the Bell measurement is given by the solution of 
Eq. [fATT]) : 



Xt{x 2 ,p 2 ) = XBm(e 2Tt x 2 ,e~ 



Pi) x 



exp j-i(l-e- Tt ) (\+n th 



After recovering the classical information, Bob performs 
on mode 2 the displacement A = g(x + ip), where g is 
the gain factor. One easily verifies that, given a state 
Pt described by the characteristic function Xt{x 2 ,p 2 ), the 
displaced state D 2 (X) p t D\(X) is described by the charac- 
teristic function 



(X2,P2) = Xt ( X2j p 2 ) e *v^A] P2 - 4 V2/m[A]^ ( A13 ) 



XD 



In order to obtain the density operator p out (or, equiv- 
alently, the characteristic function Xout{x 2 ,p 2 )), of the 
final output state of the teleportation protocol, one must 
take the average on all the possible outcomes x, p of the 
Bell measurement: 

Xout(x 2 ,p 2 ) = / dxdp V{p,x) xd(x 2 ,p 2 ) 



1 



dx dp dx" in dp i x 

ip ( V2gx 2 - x'-„ I 



exp |ii (V2gp 2 - p'[j - ip (V2gx 2 - x" n ^ } 
Xbs{x'L 0; 0,K; e-3 Tt x 2 , e^ Tt p 2 ; 0, 0; 0, 0) x 



exp 



1 



(1 



-Tt\ 



1 



)[^+nth){x 2 2 +p z 2 ) \ . (A14) 



Finally, by exploiting the relation 

(2ir)~~ 1 f dk exp{ikx} = S(x), and recalling the ex- 
pression Eq. (|A"8)) of xbs, we find that Eq. (|A14|) 
reduces to Eq. (|4]). 



APPENDIX B: HOMODYNE MEASUREMENT 
IN THE CHARACTERISTIC FUNCTION 
REPRESENTATION 

In this appendix we describe the homodyne mea- 
surement for a multimodc state in the characteris- 
tic function representation. In particular, we de- 
termine the expression for the density operator and 
for its characteristic function of the reduced state 
after post-selection, i.e. of the output state re- 
sulting from the homodyne measurements on appro- 
priately selected modes. Let us consider a three- 
mode quantum state P123, with characteristic function 
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Xi23(ai,a2,a 3 ) = Tr[p 123 D 1 (ai)D 2 (a 2 )D 3 (a 3 )} = 
Xi23(xi,Pi;x 2 ,P2;x 3 ,p 3 ), with aj = 2~ 1 / 2 {x j + ipj) 
(j = 1,2,3). Let us recall that pi 23 can be written in 
terms of the Weyl expansion: 



P123 



d 2 aid 2 a 2 d 2 a 3 XmCai, a 2 , a 3 ) x 



D 1 (-a 1 )D 2 {-a 2 )D 3 (-a 3 ) 



(Bl) 



Homodyne measurements performed on two modes, act- 
ing as two conditional measurements, reduce the three- 
mode state to a single-mode one. Let pi = p and x 2 = x 
be the results of the homodyne measurements of the 
quadratures p\ and x 2 . The homodyne measurements 
are expressed by projections on the quadrature eigen- 
states \p)\ and \x) 2 . Thus the output state p 3 subsequent 
to the conditional measurements is: 



p 3 = V- 1 (P,x)Tr 12 [\p) 1 



\X)22{X\ Pl 23 \ 



(B2) 



Here the normalization factor V(p,x) = Tr 3 [p 3 ] repre- 
sents the distribution function of the outcomes p and x. 
By using the relations 



^iW-Hxrf^) ( ( B3 ) 



2(St\Da{-a 2 )\x)2 = e~- x ^- lx ^8{x 2 ) , (B4) 
Eq. dB2j) writes 

V- X (p,x) 
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(2tt) 3 



d Xl d P2 dx 3 dp 3 e lxip - ixp2 x 



X123 (xi,0; 0, P2\ x 3 , p 3 ) D 3 (- X3 ^1 P3 j . (B5) 

Denoting by X3{x 3l p 3 ) the characteristic function of the 
state p 3 , we have 

P3 = J dx 3 dp 3 X3{X3,P3)D 3 f~ ~7^ J ■ ( B6 ) 



Comparing Eqs. (|B5|) and (|B6[) . we obtain 

X3(x 3 , P3 ) = - ^ 2 X) J d Xl dp 2 e" 1 ™ 2 x 

XmOd, 0; 0,p 2 ;x 3 ,p 3 ) • 
Moreover, it is easy to verify that 



(B7) 



V{p,x) 



(2tt) 2 

Xi23(xi,0;0,p 2 ;0,0) . 



d Xl dp 2 e lxip -' lxp2 x 



(B8) 



APPENDIX C: ANALYTICAL EXPRESSIONS 
FOR THE FIDELITIES 



using two classes of non-Gaussian entangled resources: 
The squeezed Bell-like states (J7]) and the squeezed cat- 
like states (j8]). For comparison we include also the ex- 
pressions of the fidelities associated to the Gaussian twin- 
beam and of the Buridan donkey state (jlOl The fidelity 

JF^(r, <5) for the squeezed Bell-like state (O reads: 

((l + ei~g) 2 -e ir (l-eig) 2 y x 

[A 2 - 8A(g - 1) 2 |/3| 2 + 8(g - 1) 4 |/?| 4 ] sin 2 S 

+ 2 2 T ^ - I) 2 1/3| 2 ~ A] sin 6[cos S x 
( - (1 + e^gf + e 4r (l - e* g) 2 ) + sin 5 x 
((l + e*.g) 2 + e 4 '-(l-eig) 2 )]}, (CI) 



where 



A = e- 2r - T [(l + eig) 2 + e 4r (l-eig) 2 + 

2e 2r+T (l+g 2 + 2T T , R )}, (C2) 



9 = gT, 



(C3) 



and r r ,_R is defined in Eq. ©. The fidelity T^{r, 6, 
for the squeezed cat-like state ((HJ reads: 

F ( <?Mr,S,h\) = T7 TI2 r-(cos 2 <5x 

scv ' " ,u A(l + e-H 2 sin2(5) I 

e ~i(~g-i) 2 \P\ 2 + e -h\ 2 -il(~g~m0+l3')-e r gh\] 2 x 



sin<5cos<5 ^i^~iW-n+e-gh\] 2 + C A 

+ S in 2 <5e-il«-^- er H«-^)l 2 }. (C4) 



For g = 1 and r = R = n th = 0, Eqs. JCT]) and (|C4| 
reduce to the expressions of the fidelities associated to 
the ideal protocol with unit gain, that had been originally 
determined in Ref. By letting S = in Eqs. (|Cl|) 



and (|C4|) . we recover the fidelity ^F^ B {r) associated to 



the Gaussian twin-beam resource: 



= Ji p -i(9-l) 2 |/3| 2 



(C5) 



In this appendix we report the analytical expressions 
for the fidelities of teleportation of input coherent states 



Let us notice that such a result can also be obtained by 
putting I7I = in Eq. (|C4[) . For completeness, we also 
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report the expression of the telcportation fidelity associ- 
ated to the Buridan donkey entangled resource (fTU|) : 



4lM) = |e-^-W{ 1 + £__ x 

"((l + ei^ + e^l-e^) 2 ) x 

[4(5 - 1) 2 |/3| 2 - A] + 2 cos 26 e 2r (e T g 2 - 1) x 

[A - 4(.9 - 1) 2 |/?| 2 ] - 2 S m26 (g - l) 2 (/3 2 + /T 2 ) x 

((l + eiff) a -e fr (l-ei5) J )]}. 



(C6) 
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